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We consider a single impurity with spin S embedded in a three-dimensional antiferromagnetic 
system which is close to the quantum critical point (QCP), separating magnetically ordered and 
disordered phases. Approaching the QCP from the disordered phase we study the spatial distribution 
of spin density and staggered magnetization induced by the impurity. Using two methods (self- 
consistent Born approximation and renormalization group) we found a power law decay of the spin 
density oc 1 /r 3 , and of the staggered magnetization oc 1 /r with relevant logarithmic corrections. 

We demonstrate that the local spin at the impurity site r = 0 approaches to zero at the QCP. 

We show that in the semiclassical limit of large S the problem is equivalent to the exactly solvable 
independent boson model. Our results demonstrate existence of spin-charge separation in the three 
dimensional systems in the vicinity of the QCP. 
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I. INTRODUCTION 

Quantum critical phenomena is extensively develop¬ 
ing subject in modern condensed matter physics, in both 
theoretical and experimental frontiers . 1 The most vivid 
manifestations of quantum phase transitions (QPT) arise 
in low-dimensional systems such as cuprates and iron 
pnictides. However, quantum critical behaviour is also 
found in three-dimensional (3D+time) systems. Well- 
known example of 3D compound with a magnetic quan¬ 
tum critical point (QCP) is TICUCI 3. 2 Under normal 
conditions the material is in the magnetically disordered 
phase, while pressure drives QPT to the antiferromag- 
netically ordered Neel phase. 

Quantum critical properties of a system can be sig¬ 
nificantly influenced by the presence of impurities. For 
instance, substitution of Cu atoms in the parent com¬ 
pound TICUCI 3 with low concentration of nonmagnetic 
Mg impurities creates an uncompensated spin 1/2 at the 
impurity site, which induces magnetization around the 
impurities and even lead to the formation of a long-range 
magnetic order in the macroscopic volume of the crystal . 3 
In the magnetically disordered phase the magnetization 
cloud around each impurity exponentially decays over a 
few lattice spacings from the impurity. But in the vicin¬ 
ity of the QCP the effect of impurity-induced magneti¬ 
zation can be notably enhanced. Experimental observa¬ 
tions reveal an interplay between the impurity-induced 
staggered magnetization and a quantum criticality near 
the QCP . 4 - 6 

Despite of the vast amount of theoretical work on the 
impurity-induced magnetization in quasi-ID and 2D sys¬ 
tems (see Refs . 7 - 11 and references therein), we are not 
aware of similar studies in the relation to 3D materials. In 
the present paper we consider a single impurity with spin 
S embedded in the 3D antiferromagnet (AF), which is 
close to the QCP, separating magnetically disordered and 


magnetically ordered phases. Conceptually, the problem 
is similar to the Kondo effect (see 12,13 ), since, as we show 
below, the spin cloud screens the impurity’s spin at the 
QCP (Bose-Kondo effect 14 ). We study the spatial dis¬ 
tribution of the nonlocal spin density and the staggered 
magnetization induced by the impurity. We show, that 
when approaching the QCP from the disordered phase, 
the spin density around the impurity decays as oc 1 /r 3 
with logarithmic corrections and the total spin accumu¬ 
lated in the delocalized cloud is equal to S. We also 
demonstrate that the induced staggered magnetization 
decays as oc 1 /r. 

Closely related to the problem of the impurity-induced 
spin density and the impurity-induced Neel order is a 
phenomenon of spin-charge separation (SCS). Conven¬ 
tional definition of SCS relies on the existence of two 
quasiparticles carrying spin and charge (’’spinon” and 
’’chargon”), which is the case in ID Tomonaga-Luttinger 
liquid of strongly interacting electrons. 15,16 . By contrast, 
in higher spatial dimensions there are no known sys¬ 
tems with SCS in the conventional definition. However, 
SCS exists in 2D models, such as hole-doped AF 18 - 21 . 
Furthermore, the recent research 22 reports pronounced 
SCS in the vicinity of the magnetic QCP. In the latter 
case, the precise meaning of SCS is different from SCS in 
Tomonaga-Luttinger liquid. A hole creates a spin cloud 
around the charge with radius which diverges at the QCP. 
As a result the hole’s spin becomes delocalized and spa¬ 
tially separated from the impurity’s charge pinned to the 
impurity’s site, that basically means SCS. In the present 
article we show that such SCS also occurs in 3D systems 
near the QCP. 

The paper is organized in the following way. In the 
Section II we introduce an effective field theory describ¬ 
ing 3D AF doped with a single impurity in the vicinity of 
the QCP. Considering interaction of the doped AF with 
a probe magnetic field, we introduce an operator of the 
spin density and explain how we calculate the induced 
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spin density. Here we also provide method of calcula¬ 
tion of the staggered magnetization around the impu¬ 
rity. The rest of the paper is divided into two parts, 
which correspond to the two techniques of the calcula¬ 
tions: Self-Consistent Born Approximation (SCBA) and 
Renormalization Group (RG) approach in 3 + 1 dimen¬ 
sions. Section III refers to the calculation of the spin 
density in SCBA for the most physically interesting case 
of the impurity with spin 5=1/2. In the Section IV 
we again calculate the impurity-induced nonlocal spin 
density, the local spin of the impurity and the staggered 
magnetization using RG technique. We also consider a 
semiclassical limit of the impurity with a large spin. We 
draw our conclusions in the Section V. 


II. EFFECTIVE THEORY 


An example of a 3D lattice model which incorporates 
main features of magnetic quantum criticality is pre¬ 
sented in Fig. 1. The model corresponds to the cu¬ 
bic lattice AF consisting of spins 5 = 1/2 at each site 
with weak J bonds and strong J' bonds. The system 
has the QCP driven by parameter g = J'/ J and located 
at g c = 4.013, which separates the disordered magnetic 
phase of spin dimers at g > g c from the Neel phase at 
g < <7c- 23 ’ 24 This lattice model describes various prop¬ 
erties of TICUCI 3 near the pressure-driven QCP in zero 
and non-zero magnetic field . 23 

Substitution of a 5 = 1/2 Cu 2+ ion with a spinless 
Mg 2+ creates a vacancy (hole). This is shown in Fig. 1. 
The vacancy acts as an effective impurity with the spin 
5 = 1/2. The vacancy in the lattice induces a nonlocal 
magnetization cloud around the impurity site. In the 
present paper we will calculate spatial distribution of the 
spin density and the staggered magnetization in the spin 
cloud around the impurity. 

Magnetic properties of the critical system are deter¬ 
mined by low-energy magnetic excitations. The magnetic 
excitations are magnons in the Neel phase and triplons 
in the paramagnetic phase. Hereafter we use the term 
magnons for the both types of quasiparticles. The ef¬ 
fective theory, which describes magnons in the vicinity 
of the QCP is based on the following Lagrangian, see 

25. 


e.g 


C-m = 


(d t 4>) (VA 


A 2 0 2 


,2'i2 


Qq (</> ) 
4! 
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where <j> y , (f> z ) is the magnon field, Aq oc g — 

g c is the magnon gap (squared), «o is a four-magnon 
coupling constant, dt is the time derivative, Vj = d/drt 
is the three-dimensional gradient. Hereafter we set Plank 
constant and magnon speed equal to unity H = c = 1. In 
the disordered magnetic phase Aq > 0. Near the QCP 
the magnon gap A 0 —> 0. 

The Lagrangian (1) contains quadratic terms as well 
as quartic term oc </ 4 , describing the magnon self-action. 
The magnon self-action results in the renormalization of 


the magnon gap Ao in the Lagrangian (1). From the one- 
loop RG calculations 25 it follows, that in the disordered 
phase the evolution of the renormalized gap is given by 


A 2 a A 2 
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where N = 3 in the present case of the 0(3) universal¬ 
ity class system, and C(A) is a positive constant, deter¬ 
mined by an ultraviolet scale A. Besides that, the 0 4 
term leads to renormalization of magnon quasiparticle 
residue . 25 However, the change of the residue appears 
only in the two-loop renormalization group, and there¬ 
fore is small. Hence, for our purposes we drop out the 
self-action term from the Lagrangian (1) and substitute 
the bare magnon gap to a renormalized value Ao —> A. 

The Lagrangian of a non-interacting spin-5 impurity 
reads 


C-imp = i(i>\r,t)d t ip(r,t) - (dtipi(r,t))il>(r,t)) ( 3) 

Here if> is the 25 + 1 component spinor. Hereafter we 
set the energy of the non-interacting impurity to zero. 
The Lagrangian which corresponds to the interaction be¬ 
tween the impurity and the magnon field in the disor¬ 
dered phase is 9 

Qnt =■ 4>)4>, (4) 

where A is the coupling constant, 5 = (S x , S y , S z ) are 
the operators of the impurity spin acting in the (25 + 1)- 
dimensional Hilbert space. 



Figure 1: An example of a lattice model for 3D AF with 0(3) 
QCP. Spins S = 1/2 located at each site. Thin lines denote 
weak J bonds and thick lines denote strong J' bonds. A 
quantum phase transition between the Neel and the dimerized 
paramagnetic phases occurs at (J'/J) c = 4.0 13. 23,24 The big 
blue sphere represents an impurity (hole) introduced into the 
lattice. 


The interaction of the impurity with magnons leads to 
appearance of the nonlocal part of spin density s{r). In 
order to find s(r) we use the Lagrangian of interaction 
of the system with external magnetic field (see, e.g. 26 ) 


C B = - ■ [B x 0 ]) + 


[B x ft 


(5) 
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We have set here hbQ = 1- Note, that Eq. (5) is still 
valid, if the magnetic field B is non-uniform. In the 
contrast to the majority of previous works, where B is 
considered uniform, the non-uniformity of the probe mag¬ 
netic field is crucial for the present paper. Linear in B 
terms in the Lagrangian (5) provide following expression 
for spin density: 

s(r) = ([<t> x dt<f>\ + h.c.) + ^Sip) = 

s n i(r) + S imp S{r). (6) 


The brackets (• • •) denote an averaging over the ground 
state of the system. The term 1/2([</> x dt<fi\ +h.c.) in Eq. 
(6) is the nonlocal part of the spin density s n i (r ), induced 
by the impurity. The subscript ”nl” stands hereafter for 
’’nonlocal”. The term (i/.d Si/j) in Eq. (6) corresponds to 
the local spin Si mp at the impurity’s site. 

In addition to the spin density we will consider the 
staggered magnetization, induced by the impurity. Writ¬ 
ing down an Euler-Lagrange equation for the magnon 
field cj) from the action f dtd 3 r {Cm + £mt} and tak¬ 
ing expectation value of the result we obtain Yukawa-like 
form of the staggered magnetization 


(4>{r)) = -A- 


„ —A r 


Airr 


( 7 ) 


At the QCP the exponent in Eq. (7) is close to unity and 
(</>(r)} oc 1/r. Therefore, in order to find corresponding 
prefactor, we only need to calculate the local spin at the 
impurity site Si mp . 

To find local as well as nonlocal components of the 
impurity-induced spin density we will proceed with the 
following procedure. We calculate the shift eg of the 
ground state energy, corresponding to the probe mag¬ 
netic held B(r') = B8(r' — r). The energy shift of the 
system reads = B ■ s(r), therefore 


s(r) 


de B 

dB 


( 8 ) 


The spin density s(r) = es(r) and the staggered mag¬ 
netization (0(r)) = e{(j>{r)) are directed along the impu¬ 
rity spin Si mp = e Sim P (e is a unit vector), and due to 
spatial isotropy of the system depend only on r = |r|. 
The ground state energy e g of the system is the position 
of a singularity of the retarded impurity’s Green’s func¬ 
tion Gg(e) and can be found from the Dyson’s equation 

Gs\e)=e-t(e)-B^(e,r)=0, (9) 


where £(e) is the self-energy of the impurity at zero mag¬ 
netic Held, T /i (e, r) is the vertex function, corresponding 
to the interaction of the system with the probe magnetic 
Held. Note, that in Eq. (9) we need to keep only linear 
in terms. From the rotational symmetry properties 
the only possible ’’kinematic” structure of the vertex 

f^(e,r)=T(e,r)>/5. (10) 


(b) 



Figure 2: Example of diagrams for a coupling between a probe 
magnetic field B and (a) impurity spin, (b) nonlocal spin 
density. Solid line corresponds to a bare (A = 0) impurity 
Green’s function, waivy line represents magnon, dashed line 
represents probe magnetic field B. The cross on the magnon 
line corresponds to the magnon-B vertex, provided by the 
term B[(f> x cj>\ in the Lagrangian (5). 


The vertex function can be split in local and nonlocal 
parts 


f^( e ,r) 


rf mp (e), r = 0, 

r{Ji(e,r), r>0. 


( 11 ) 


Calculating the shift es of the position of the singularity 
in the Green’s function G B {e) due to the probe magnetic 
field and using the formula (8) we find the local and non¬ 
local components of spin density s(r). 

Below we will calculate spin density using two ap¬ 
proaches: Self-Consistent Born Approximation (for S = 
1/2) and Renormalization Group (for arbitrary S). 


III. SELF-CONSISTENT BORN 
APPROXIMATION (S=l/2) 


Standard approach for calculation of a single-fermion 
Green’s function is 1/A expansion for O(N) group, where 
N = 3 is the number of magnon components. Summa¬ 
tion of leading terms in the 1/N expansion results in Self- 
Consistent Born Approximation (SCBA), see Fig. 3. We 
will apply SCBA to the case of S = 1/2 impurity only. 
As it will be demonstrated in Section IV C, for S > 1/2 
corrections to impurity-magnon vertex, which are disre¬ 
garded in SCBA, become relevant. Therefore, in this case 
SCBA fails and application of RG technique is necessary. 


A. Impurity Green’s function at zero magnetic 
field 

To consider interaction of the system with the probe 
magnetic field, we first calculate Green’s function of the 
impurity at zero magnetic field. The Green’s function of 
the impurity at B = 0 is proportional to the identity ma¬ 
trix in the spin space G(e) = G(e). The Dyson’s equation 
for the Green’s function is graphically represented in Fig. 
3. The analytical form of the equation reads: 

G(e) = 


e — £(e) + *0 ’ 


( 12 ) 
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where the impurity self-energy is given by following ex¬ 
pression 

t(e) = X 2 J = 

3 £ M lG(e - uj q ). (13) 

Q 

Here uj q = -\/A 2 + q 2 is the magnon dispersion, M q = 
X/■ s J2uj q is the matrix element corresponding to emis¬ 
sion of magnon with momentum q by the impurity, 
q) = S I1V / (w 2 — w 2 + iO) is the magnon propaga¬ 
tor. We expressed spin 1/2 operators via Pauli matrices 
S M = cr M /2. Combinatorial factor 3 in Eq. (13) comes 
from summation over intermediate polarization state of 
magnon. 



Figure 3: Dyson’s equation in SCBA. Double line is impurity’s 
Green’s function. 


The sum over momentum q in Eq. (13) diverges at 
large |q|, therefore we have to introduce ultraviolet cutoff 
A. The parameter A depends on particular realisation of 
system and can be estimated as inverse lattice spacing in 
the host AF. 

Solution to Dyson’s equation (12) near the QCP (A —> 
0) has the following form 


G ” 1(e) - (£ “ 60 + *° } \/ 1 + ^ ln ( eo + A A c -*o) 

(14) 

in the vicinity of the singularity point eo ~ — 3AA 2 /47r 2 . 
Formula (14) is obtained with logarithmic accuracy, i.e. 

assuming that In 1. 

The Green’s function (14) has nontrivial analytic struc¬ 
ture. At finite magnon gaps the quasiparticle pole at 
e = 6q is separated by the gap A from the incoherent 
part of the Green’s function. At the QCP, when A —> 0 
the pole and the branching point singularity are merg¬ 
ing. The quasiparticle residue of the impurity Green’s 
function G(e) vanishes in the vicinity of the QCP 


Z = 


SEMy 1 

de ) 


\A + §^ ln (i) 


o. 


A—>0 


(15) 

Vanishing quasiparticle residue is the first signal of delo¬ 
calization of the impurity-induced spin cloud and there¬ 
fore indication of SCS. 22 

Typical value of the impurity-magnon coupling con¬ 
stant A can be estimated on the basis of the lattice model, 


shown in Fig. 1. Lattice calculations 27 result in the value 
of the effective coupling constant n = 3A 2 /27t 2 ~ 0.5, ap¬ 
pearing in front of the logarithm in formula (14). There¬ 
fore, the logarithmic corrections are significant in the 
vicinity of the QCP. 

Analytical result (14) for the impurity Green’s func¬ 
tion can be compared with direct numerical solution of 
the Dyson’s equation (12), corresponding plots for spec¬ 
tral functions — l/7rIm{G(e)} are plotted in Fig. 4. An 
artificial broadening iO -A «2.5 x 10~ 3 A is introduced in 
the numerical procedure and in analytical formula (14). 
We see an excellent agreement between the numerical and 
the analytical results. 

Let us make a comment about the validity domain of 
SCBA for the results in the present section, and all fol¬ 
lowing results, which will be derived in Sections IIIB, 
IIIC. Formally, SCBA relies only on 1/N expansion of 
O(N) group, independently on the value of the coupling 
constant A. SCBA is applicable for arbitrary A, in con¬ 
trast to RG method, which works only for A < 1. We 
will return to this discussion later, in the Section IV A. 



Figure 4: Spectral function of S = 1/2 impurity obtained in 
SCBA. Panel (a) corresponds to the QCP (A = 0); panel (b) 
corresponds to magnon gap A = 0.05A. Effective coupling 
constant is set to k = 0.6. Solid black line corresponds to 
Green’s function, calculated numerically, red dashed line cor¬ 
responds to analytical formula (14). Note, that on the panel 
(a) the position of pole and branching point are merging. 


B. Calculation of nonlocal spin density s n i(r) 


To evaluate nonlocal spin density induced by the im¬ 
purity at the distances r > 0, we substitute the Green’s 
function pole position e g = eo + cb into Eq. (9), expand 
it in £b up to the first order and use Eq. (8). The result 
reads 


s n i(r) = ZT n i(e 0 ,r ). (16) 


Leading in coupling constant A contribution to the ver¬ 
tex Ki is represented by the Feynman diagram shown in 
Fig. 2, (b). The analytical expression for the diagram is 
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following 


r^(e,r) =T nl (e,r)a lx = f ^ (Acr^) G 0 (e - w) 


(A aP)D„ u (tu,k) 


2zCJ £ fii/a 6 


q.k 
i(q—k)r 


Dp a (uj,q) 


(17) 


where Go(e) = l/(e + iO) is the bare retarded Green’s 
function of a non-interacting impurity. Expression in 
square brackets corresponds to magnon - probe magnetic 
field vertex, which we show in Fig. 2, (b) as a circle with 
a cross inside. 

SCBA equation for vertex r^(e,r) is graphically rep¬ 
resented in Fig. 5. Analytical form of the equation for 
the vertex is 

T n i(e,r) =r^ ) (e,r) - 
Y,M 2 q G 2 (e-u; q )r nl (e-Lu q ,r), (18) 

Q 


where r^(e, r) corresponds to the first term in the rhs 
of the diagrammatic equation in Fig. 5 and reads 


r<°>(e,r) = 2A 2 £\ 


— k)r - Vq) - G(e - UJk ) 


,i(q-k)r 


q.k 


, ,2 , ,2 


.(19) 


To obtain expressions (18) and (19) we performed inte¬ 
gration over u in the rhs of the original SCBA equation, 
shown in Fig. 5. Factor (—1) in Eq. (18) comes from 
algebraic identity for Pauli matrices cr^cAcr^ = — cr v . 
Fformula (19) follows from (17), where the bare impu¬ 
rity Green’s function is changed to the ” dressed” Green’s 
function G(e), shown in Fig. 3. 



Figure 5: Diagrammatic equation for ’’nonlocal” vertex func¬ 
tion f£,. 

To find s n i(r) we solve numerically Eq. (18) for the 
vertex T n i(e,r) and substitute the result together with 
the quasiparticle residue Z, obtained from numerical so¬ 
lution of Dyson’s equation (12), to Eq. (16). Solution to 
Eq. (18) has been found iteratively, starting iterations 
from the T ni (e,r) =r^ } (e,r). The results of the calcula¬ 
tion of the spin density s n i(r ) for different values of the 
magnon gap A and coupling constant A are presented in 
Fig. 6. 

For the purpose of computational efficiency we used 
spherical cutoff |q|,|fe| < A in integrals in Eq. (18) 
and (19), instead of integrating over a cubic Brillouin 
zone. This cutoff scheme resulted in appearance of sig¬ 
nificant r-oscillations in the induced spin-density s n i(r), 



Figure 6 : Spin density s n i(r) (multiplied by r 3 ) induced 
by S' = 1/2 impurity as a function of dimensionless dis¬ 
tance y = Ar calculated in SCBA. Points represent nu¬ 
merical results for different values of magnon gap and effec¬ 
tive coupling constant k = 3A 2 /27r 2 . Squares correspond to 
A = 6.25 x 10 _ 3 A, triangles correspond to A = 1.25 x 10 _ 2 A; 
filled markers represent k = 0.3, open markers correspond to 
k = 0.6. Solid line is analytical approximation (20). 


where the period of oscillations is r ~ 1/A and the ampli¬ 
tude of the oscillations decaying with increasing r. It is 
clear that these oscillations are by-products of the rigid 
spherical cutoff and will be notably suppressed, if one 
performs proper 3D-integration over the cubic Brillouin 
zone. Hence, in Fig. (6) we plot numerical data for the 
spin-density, averaged over the period of the oscillations. 


Our numerical calculations show that the starting ap¬ 
proximation r ( n ° } (e, r ) for vertex function and the solution 
T n j(e, r ) °f the SCBA equation (18) are very close to each 
other. Therefore, to obtain analytical approximation for 
the nonlocal spin density we substitute (19) in Eq. (16) 
and use formula (14) for the impurity’s Green’s function, 
the result reads 

, , A 2 A iG(2Ar) 

Wl + f^n^l + f^lnAr ^ 

( 20 ) 

Here K i(x) is Macdonald function of the first kind. At 
distances 1/A < r < 1/A using expansion of the Mac¬ 
donald function K\(x) —> \/x at x —^ 0, we obtain power- 
law asymptotics with logarithmic corrections for the spin 
density: 


s n i(r) 


A 2 


87 r 3 r 3 


\/(f + f^ hl s) I 1 + H? lnAr ) 


■ ( 21 ) 


At large distances r > 1/A the spin density (20) is ex¬ 
ponentially suppressed: s n i(r) oc e~ 2Ar /r 5//2 . In Fig. (6) 
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field. The result for the impurity spin reads 




imp — ^ rimp(^) 


(23) 


Diagrammatic equation for the vertex function ff mp (e) = 
ri m j,(e)cr M in SCBA has the graphical representation, 
shown in Fig. 8. Corresponding analytical form of the 



Figure 8: Diagrammatic equation for ’’local” vertex function 
f ? . Cross represents bare vertex = 1. 


Figure 7: Integral nonlocal spin S„i and local spin of the im- 
putity Simp as a function of the magnon gap A in SCBA. 
Full and open markers correspond to values of effective cou¬ 
pling constant k = 0.6 and k = 0.3. Circles correspond to 
Sni, squares represent Simp and triangles show the net spin 
Sni T Simp- Solid lines are theoretical predictions for the local 
impurity spin Simp , given by Eq. (25). Dotted lines are visual 
guides for S n i. Red dashed line corresponds to the net spin 
equal to 1/2. 


solid lines correspond to the analytical result given by 
Eq. (20). One can see excellent agreement between the 
analytical and the numerical results. 

The net spin of the system, which is given by the sum 
of local impurity spin and spin of nonlocal cloud is con¬ 
served and must be equal to S = 1/2. The integral spin, 
corresponding to the nonlocal spin density 

S n i = J d 3 rs n i(r) (22) 

is plotted in Fig. 7 versus A/A. We use the numerical 
results for s n i(r), shown in Fig. 6, in order to obtain S n i. 
One can see that the nonlocal spin S n i logarithmically 
increases with decreasing A and tends to S n i = 1/2 at 
the critical point. Therefore the rest of the spin should be 
attributed to the impurity spin Si mp = 1/2 — S n i, which 
vanishes at the QCP. We check this statement in Section 
IIIC, calculating local spin of the impurity. 


C. Local spin of the impurity and staggered 
magnetization 


To calculate impurity spin, localized at r = 0, we use 
similar approach to that we used in the previous section. 
We introduce local magnetic field B(r) = B5(r ) and 
calculate energy shift of the impurity due to the magnetic 


equation represented in Fig. 8 is 

r imp (e) = 1 -^2 M q G 2 (e - oj q )T lmp (e - oj q ) . (24) 

q 

Solving equation (24) and substituting the solution in to 
Eq. (23) we obtain with logarithmic accuracy spin of the 
impurity 


Simp — 


2(i + !Si »!) 2,s ' 


(25) 


We also calculate residue spin of the impurity numeri¬ 
cally, solving iteratively Eg. (24). Both analytical and 
numerical results for Si mp at different values of the pa¬ 
rameters A and A are plotted in Fig. 7. We see good 
agreement between the analytical and the numerical re¬ 
sults. From Fig. 7 we can notice that the impurity spin 
logarithmically tends to zero when we approach to the 
QCP, A —> 0. In Fig. 7 we also show the net spin 
Simp + S n i for different values of A and A. 

Results of our calculations, presented in Eq. (21), Eq. 
(25) and also in Figs. 6, 7 show that at the QCP the lo¬ 
cal spin is approaching to zero and the spin of the system 
is accumulated in the nonlocal spin cloud. This delocal¬ 
ized spin cloud around impurity has size proportional to 
inverse magnon gap r — 1/A, and therefore significant 
part of impurity spin is separated from charge, localized 
at r = 0. We will return to this discussion again in the 
Section IV B. 

The net spin of the system equals to 1/2. This is an 
exact statement and can be demonstrated at the dia¬ 
grammatic level. One can trace mutual cancellations of 
corrections to the impurity spin and integral spin of non¬ 
local cloud in every order in A. Corrections to the impu¬ 
rity spin Si mp are cancelled by corrections to integral spin 
S n i ■ The numerical results for the net spin of the system 
presented in Fig. 7, are consistent with the conservation 
of spin. 

Using Eq. (7) and Eq. (25) we obtain following ex¬ 
pression for the staggered magnetization induced by the 
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spin 1/2 impurity 


-Ar 


(</(r)) = -A- 


47rr /'i i 3A 2 A\ 2 / 3 


(1 + 


2^ hl A, 


(26) 


Away from the QCP the staggered magnetization, in¬ 
duced by the impurity is exponentially small. In the 
vicinity of the QCP prefactor in Eq. (26) becomes log¬ 
arithmically suppressed, however the staggered magneti¬ 
zation decays only as (<j>(r)} oc 1 jr. 


IV. RENORMALIZATION GROUP APPROACH 
IN 3+1 D 

In this section we calculate nonlocal and local compo¬ 
nents of the spin density using RG technique in 3 + 1 
dimensions. In the RG approach the coupling constant A 
becomes dependent on the energy scale. Since 3 + 1 D is 
the upper critical dimension, the evolution of the running 
coupling constant is logarithmic. It leads to logarithmic 
corrections to s n i and (</>(r)), similar to results (21) and 
(26) obtained in SCBA. We derive our results for the case 
of an arbitrary spin S of the impurity in Sections IV A 
and IV B, and then we analyze the limit of a large spin 
S in Section IV C. 

In RG technique we consider evolution of the coupling 
constant A, quasiparticle residue Z, spin density and 
staggered magnetization with the energy scale /i, start¬ 
ing evolution from the ultraviolet scale A and finishing at 
the infrared scale A. The scale /i here has the meaning 
of the characteristic energy transfer from magnons to the 
impurity. At ultraviolet scale A we set parameters of the 
theory to bare values, in our calculations A plays a role 
of a renormalization point. Observables in the vicinity 
of the QCP are calculated as the result of RG evolution 
from the ultraviolet scale A to the infrared scale /r = A. 


A. Evolution of coupling constant and 
quasiparticle residue 

First, we calculate evolution of the coupling constant 
A(/r). One-loop correction to the coupling constant is 
represented by the sum of diagrams, shown in Fig. 9. 



Figure 9: One-loop corrections to the impurity-magnon cou¬ 
pling constant A. 

Note, that in RG approach the correction to the cou¬ 
pling constant includes the vertex correction [Fig. 9, (b)], 
and also the self-energy correction [Fig. 9, (c)]. This is 
different from SCBA, in which we disregard diagram (b). 


Contribution dA (b ' to the coupling constant correction 
is given by the diagram (b) in Fig. 9 and reads 


QIS Q/J, QU 

S'* 6 A (b) = A 3 ——x 


S 2 


iduj r 
27r 


G 0 (n - u')G 0 (p, - w - u')D(u', k) 


SU (S(S + 1)-1) A 3 h A 


S 2 


47T 2 


(27) 


After cancelling out the factor S M from the both sides of 
Eq. (27) we obtain SX^ b \ The second contribution SX^ C> , 
which comes from the diagram (c) in Fig. 9 reads 


<5A (c) 


G 0 (/i- uj) x 

/ ~ w “ u') D (u)' ,fc) 



(28) 


and contains linear in A and logarithmic in A terms. The 
linear term corresponds to the shift of the position of 
quasiparticle pole eo in the impurity’s Green’s function 
and therefore is irrelevant for our purposes. The loga¬ 
rithmic term in <5AG) reads 

{a "’- > -;£( 1 + !) 1 V (29 > 

Total correction to the coupling constant A is 

SX = SX^ + 5X^ = --A— In -. (30) 

4S 2 tt 2 fi y ’ 

Note, that for S' = 1/2 the vertex correction (27) is sup¬ 
pressed by the factor 1/N = 1/3, comparing to SX^. 
This suppression corresponds to standard 1 /N expansion 
of O(N) group. However, at large S the 1/N suppression 
of SX <h '> is compensated by S, and hence S A® and <5A( C ) 
to a large extent compensate each other, SX <h ^ ss — SX^. 
Thus, at large S the vertex correction becomes signifi¬ 
cant and can not be disregarded. This is the reason why 
SCBA fails in the case of large impurity spin. 

In the paradigm of RG, evolution of physical parame¬ 
ters on some energy scale ^ is determined by the value 
of X(n) on the same scale. Hence, Eq. (30) results in the 
following Gellman - Low equation 

dXjfJ,) _ A 3 (/z) 

din/i 4 j5' 2 7t 2 1 j 

Solution to Eq. (31) with the initial condition A (A) = A 
is 


AM = : \ 2 A ■ (32) 

V 1 + 25 2 7r 2 II 

Note, that the running coupling constant (32) vanishes 
in infrared limit: A(/i) —>Oat/i~A—>0. The RG 













scale /r is bounded from below by the value of magnon 
gap fi > A. 

In order to find the quasiparticle residue of the impu¬ 
rity’s Green’s function we consider one-loop correction to 
the impurity’s self-energy. Logarithmic part of this cor¬ 
rection was already calculated as a part of the diagram 
(c) in Fig. 9. Corresponding equation for evolution of 
Z(jj,) reads 


din Z(n) = / 1'\ A 2 (^) 

d In fj, \ S) 47t 2 


(33) 


Solution to Eq. (35) with initial condition S/ mp (A) = S 
is following 


Simp (/0 


s 



A 2 

2S 2 7r 2 



(36) 


and proportional to solution (32) for running coupling 
A (/Li). The local spin at the impurity site is equal to 
Simpifi — A) and approaches to zero at the QCP. Using 
the result (36) and relation (7) we obtain distribution of 
staggered magnetization around impurity 


Solution to Eq. (33) with initial condition Z(A) = 1 
reads 




1 



A 2 

2S 2 7r 2 



S(S+l)/2 


^A(/r)^ S(S+1) 


(34) 

The quasiparticle residue Z(/i) vanishes, while approach¬ 
ing to the QCP: /x ~ A —> 0. 

Note, that RG approach, being used in the current 
Section, is valid if the effective coupling constant k = 
X 2 /2S 2 n 2 < 1, since we perform perturbative expan¬ 
sion, such as in Eq. (30). However, the proper expan¬ 
sion parameter in the vicinity of the QCP is not k , but 
rein A/A. RG method (in single-loop approximation) al¬ 
lows to sum up (leading) logarithmic corrections of the 
following kind k m ^2 n k n In" ( A/ A) . Therefore, the re¬ 
sults obtained within one-loop RG in the Section IV are 
valid when k < 1 and K 2 ln(A/A) < 1, but the product 
k1ii(A/A) can have an arbitrary value. 


B. Impurity spin and nonlocal spin density 

Now we consider RG evolution of the impurity spin 
Si mp and spin density distribution s(r) with renormal¬ 
ization scale /x. As in Section III we calculate S, mp and 
s(r), considering interaction of the system with probe 
magnetic field B{r). 

We start from calculation of corrections to S lmp due 
to interaction of the impurity with magnons. One-loop 
corrections to Si m p are shown in Fig. 10. 

I + I 

' ' (°) ' ( & ) ' (c) 

Figure 10: One-loop corrections to "local” spin Simp. 



(<^(r)) = 


A e~ Ar 

47IT R ' A2 n, A 

V 1 + 2S 2 7r 2 ln A 


(37) 


Now we calculate the nonlocal spin density s n i(r). In 
RG technique it is more natural to use momentum rep¬ 
resentation for the spin density, therefore we write evolu¬ 
tion equation for the Fourier component s n i ( q). Leading 
in A 2 contribution to s n i(q) is provided by one-loop dia¬ 
gram, shown in Fig. 2, (b). Evaluation of this diagram 
with logarithmic precision leads to 


s ( o) 

^nl 


(9) 


^ In A M > q, 

4^1nf> ^<9- 


(38) 


Fourier transform of the second line of Eq. (38) gives 
the spin density s^(r) = X 2 /16Sn 3 r 3 at the distances 
1/A < r < 1/A. In analogy with the result (21), ob¬ 
tained in SCBA, in RG calculations, we should expect 
logarithmic corrections to oc 1/r 3 distribution. Note, 
that the logarithmic corrections are important, because 
they provide proper normalization condition of the inte¬ 
gral nonlocal spin f d 3 rs n i(r ) —► S at the QCP. Volume 
integral of the spin density oc 1/r 3 is logarithmically di¬ 
vergent oc In A/A if we disregard the log corrections. 

In order to account for RG evolution of the spin- 
density, we evaluate single-loop corrections to the leading 
diagram presented in Fig. 2, (b). Diagrams (b) and (c) in 
Fig. 11 represent these corrections, which are similar to 
corresponding diagrams in Figs. 9 and 10. RG evolution 



(a) ' (6) ' (c) 


0 

I 


Note, that the diagrams in Fig. 10 are analogous to 
the diagrams in Fig. 9 for corrections to the coupling 
constant A. The only difference is that the impurity- 
magnetic field coupling Si mp • B is proportional to the 
impurity spin Si mp . Hence, RG evolution equation reads 


dSimpiR) _ A (/x) G ^ 

dinii ~ AS 2 n 2 impW 


( 35 ) 


Figure 11: One-loop corrections to nonlocal spin density 

of nonlocal spin density distribution reads 

ds n i(q, /x) _ j ±^s n i(q,[i) - M»9, /one 

\ js$s n i(q,IJ.), /!<<?• 
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Note, that equation for evolution of the spin density 
with p, is different in two domains (i > 9 and /Li <gC < 7 , 
which is due to the fact, that one-loop expression (38) 
for spin density has various form in the both domains. 
We solve Eq. (39) separately in the two domains and 
match the solutions at fi ~ q. As an initial condition 
for the evolution equation (39) we set s n i( A) = 0. We 
obtain following result for the spin density at infrared 
scale fx — A: 


Snl(q) 


S 


1 + 


A 2 

25^7r 


■ In 


A 



( 4 °) 

Condition of the net spin conservation in the momentum 
representation has the form s n i (q) | 0 +5j mp = 5. Using 
expressions (36) and (40) it is easy to check the net spin 
conservation, having in mind that the low bound for the 
momentum q in our formulas is q ~ A. 

Calculating Fourier transform of Eq. (40) we obtain 
spatial distribution of the induced spin density 


'-’ru ) — j - j - 

165 7 r 3 r 3 - v /l+ 2 ^rln^yi+ 2 ^ 1 nAr 

(41) 

at the distances 1/A < r < 1/A. Using (41) and 
(36) one can verify conservation of the net spin in r- 
representation: f d 3 rs n i(r) + Si mp = S. Integration of 
nonlocal spin density should be performed in the range 
of the distances 1/A < r < 1/A which is defined by the 
infrared and ultraviolet cutoffs for our theory. 

Note, that at the QCP the main contribution to the 
nonlocal spin f d 3 r s n i(r ) comes from large distances r < 
1/A —> 00 . Indeed, the integral 


C. Semiclassical limit: impurity with large spin 

From theoretical point of view it is interesting to con¬ 
sider the semiclassical case of large spin of the impurity. 
Taking formal limit 5—1 00 in Eqs. (36), (37) and (41) 
we obtain 

Sim P = S , s n i{r) = 0, (0(r))=--. (43) 

47rr 

We see from Eq. (43), that in the semiclassical limit 
there is no nonlocal spin density around impurity and 
local spin S is ’’unscreened” in this case. Therefore, there 
is no spin-charge separation in the semiclassical limit. 

Note, that local impurity spin, nonlocal spin density 
and staggered magnetization in the semiclassical limit are 
provided just by tree-level approximation. The reason is 
that quantum fluctuations of the of impurity spin are 
suppressed at large S. Indeed, let us consider the case of 
impurity in the state with maximal projection of spin on 
the quantization axis z: \ 5, S z = S). Interaction of the 
impurity with magnon either leaves projection S z to be 
unchanged or changes it by unity, A S z = — 1. The action 
of operator 5~ on the state |5, 5) provides eigenvalue 
S. On the other hand, matrix element of the lowering 
operator 5_ between states 15, ^—1) and | S, S) is equal to 
a/ 25- Therefore, the processes with change of projection 
of the impurity spin are suppressed in the limit of large 
5. 

In the semiclassical limit only z-component of the op¬ 
erator of impurity spin is relevant, as a result the La- 
grangian in Eq. (4), corresponding to interaction of the 
impurity with magnons takes the form 

Cint = (44) 


/ 

l/A <r<R 


d 3 r s n i(r) = S 


\A + 2 l n Ai? - 1 
\J l + 23^ ln S 


(42) 


logarithmically grows as a function of the upper integra¬ 
tion limit R, which means that major part of spin in the 
nonlocal cloud is accumulated at the distances of the or¬ 
der of R ~ 1 /A. At the same time, the local spin of the 
impurity Si mp vanishes at the QCP, see Eq. (36). There¬ 
fore, we conclude that at the QCP the impurity’s spin is 
spatially separated from the impurity charge. 

The results (34), (37) and (41) obtained in RG tech¬ 
nique are similar to corresponding answers (15), (26) and 
(21), obtained in SCBA. For the spin 5=1/2 the differ¬ 
ence is in the numerical factors in front of the logarithms: 
3X 2 /2n 2 in SCBA, comparing to 2A 2 /n 2 in RG. For Z 
and (<j>{r)) the powers of logarithms are also insignifi¬ 
cantly changed: 1/2 —1 3/8 and 2/3 —>- 1/2, respectively. 
The reason for this minor changes is due to the 1/A ver¬ 
tex correction, which is accounted in RG approach (see 
diagram (b) in Fig. 9), and is disregarded in SCBA. The 
RG results are more accurate, then SCBA results. How¬ 
ever, the expansion of the RG results and SCBA results 
coincide up to the single-loop order (first order in A 2 ). 


Thus, the problem of classical impurity ’’dressed” with 
z-polarized magnons is equivalent to the problem of im¬ 
purity interacting with scalar bosonic field <f> z . 

The problem of interaction between impurity and 
scalar boson field is known as independent boson model, 
this model is exactly solvable . 28 The exact solution agrees 
with Eq. (43). 

Retarded Green’s function of the impurity in time rep¬ 
resentation at t > 0 reads 28 


G(t) = —iex p 


ite 0 — A 2 ^ 

Q 


(1 - e~ iUqt ) 
2 uj 3 


(45) 


where eo = —A 2 l/2w 2 . Performing Fourier transfor¬ 

mation of the impurity Green’s function (45), and calcu¬ 
lating quasiparticle residue at the Green’s function pole 
e = eo, we obtain 



In the limit 5 —> 00 the RG result (34) is consistent with 
Eq. (46). 
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V. DISCUSSION OF RESULTS AND 
CONCLUSION 

In the present paper we have considered a single im¬ 
purity with spin S embedded into 3D AF system, which 
is close to the 0(3) quantum critical point (QCP), sep¬ 
arating paramagnetic and Neel magnetic phases. The 
impurity spin induces the usual magnetization and the 
staggered magnetization clouds around position of the 
impurity. Using the effective Lagrangian method and ap¬ 
proaching the QCP from the disordered phase, we have 
calculated spatial distributions of the spin density s(r) 
(magnetization) and the staggered magnetization (< j>(r )) 
in the cloud. For calculations we use two different meth¬ 
ods, Self Consistent Born Approximation (SCBA) and 
Renormalization Group (RG). SCBA is justified by the 
small parameter 1/N where N = 3 for the 0(3) group, 
while RG is justified by the small coupling constant. We 
show that for S' = 1/2 results of both methods are consis¬ 
tent within expected accuracy 1/N. However, at larger 
values of the impurity spin the SCBA method is not valid 
because the small parameter 1 /N is compensated by the 
large spin. Therefore, for S > 1 only RG results are valid. 

The impurity quasiparticle residue vanishes at the 
QCP, see Eq. (34). This is the first indication, that 
the impurity spin is fully transferred to the magnon 
cloud. The effect of screening of impurity’s spin by 
spin-one magnetic fluctuations is a Kondo effect in a 
bosonic sector. 14 The spin density has a local component 
SimpS(r), which is localized at the site of the impurity, as 
well as a spatially distributed nonlocal part s n i(r). Be¬ 


cause of the vanishing residue the average impurity spin 
Simp logarithmically vanishes at the QCP, see Eq. (34). 
Of course, the total spin S is conserved and it is trans¬ 
ferred into the nonlocal spin cloud. The nonlocal spin 
density at r < 1/A, where A is the magnon gap, decays 
as s n i(r) oc 1/r 3 with proper logarithmic corrections, see 
Eq. (41). Obviously, at r > 1/A the spin density decays 
exponentially. 

Integral spin in the nonlocal spin cloud is mainly accu¬ 
mulated at large distances r ~ 1/A, see Eq. (42). There¬ 
fore the spin is spatially separated from the impurity and 
at A —> 0 the separation scale becomes infinite. In this 
sense our results demonstrate the spin-charge separation 
in 3D magnetic systems at the QCP. 

Interestingly, the cloud of the staggered magnetization 
at r < 1/A decays only as the first power of distance, see 
Eq. (37). This is why a tiny concentration of impurities 
can significantly influence the critical behaviour of the 
system. 

Finally, we have analyzed the semiclassical limit of a 
very large impurity spin, S 3> 1. In this limit the quan¬ 
tum spin-flip transitions become negligible and the spin 
impurity problem is reduced to an exactly solvable text¬ 
book example. 28 
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